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Abstract. The difference in stellar structure above and below spectral type ~M4 is 
expected to be a very important one, connected directly or indirectly to a variety of 
observational phenomena in cool stars — such as rotation, activity, magnetic field gen- 
eration and topology, timescales for evolution of these, and even the basic mass-radius 
relationship. In this Cool Stars XVI Splinter Session, we aimed to use the M4 transition 
as an opportunity for discussion about the interiors of low-mass stars and the mecha- 
nisms which determine their fundamental properties. By the conclusion of the session, 
several key points were elucidated. Although M dwarfs exhibit significant changes 
across the fully convective boundary, this "M4 transition" is not observationally sharp 
or discrete. Instead, the properties of M dwarfs (radius, effective temperature, rotation, 
activity lifetime, magnetic field strength and topology) show smooth changes across 
M3-M6 spectral types. In addition, a wide range of stellar masses share similar spec- 
tral types around the fully convective transition. There appears to be a second transition 
at M6-M8 spectral types, below which there exists a clear dichotomy of magnetic field 
topologies. Finally, we used the information and ideas presented in the session to con- 
struct a framework for how the structure of an M dwarf star, born with specific mass 
and chemical composition, responds to the presence of its magnetic field, itself driven 
by a feedback process that links the star's rotation, interior structure, and field topology. 
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1. Introduction 

The interior structure of a star is a function of its mass. Above a threshold mass of ~0.3 
Mq, the interior is expected to resemble that of the Sun, having a radiative zone and a 
convective envelope. In contrast, below the threshold mass the convection zone extends 
all the way to the core. 

Despite not having a strong radiative/convective interface (like that found in the 
Sun), mid- to late-type M dwarfs (spectral type of M4 and greater) are observed to 
have strong magnetic fields and subsequent surface heating that results in observed 
magnetic activity (traced by strong line or continuum emission). Recent theoretical 
modeling of stars with convective envelopes that extend to the core (and a lack of a 
strong radiative/conv ective interface) were able to produce strong, long lived, large 
scale magnetic fields ( Brownin^l2008h . Observational results also suggest that in mid- 



to late-type M dwarfs, the fraction of the magnetic fi eld in large scale components is 
larger than in early type M dwarfs (Donati et al.. 20081 : iReiners & Basrill2008.) . 



These are just some examples of the phenomena and questions that play out around 
the M4 spectral type. The purpose of this Cool Stars XVI Splinter Session was to 
use this specific transition point — the M4 transition — as an opportunity for focused 
discussion and "out of the box" collective thinking about the interiors of low-mass stars, 
what we know, what we don't know, what we wish we understood better. We intended 
this to be a rich, fun session to bring together observational evidence from multiple 
directions with the goal of disentangling the effects of multiple physical phenomena 
and leading toward a clearer set of guiding questions for future work. 



2. The M4 transition 

Neill Reid began by reminding us of some of the basic properties of M dwarfs and of 
the basic predictions of theoretical stellar evolution models for M dwarf interior struc- 
ture. Most theoretical stellar evolution models predict the transition into the fully con- 
vective regime to occur at a stellar mass of ~0.3 Mq, but differ in detail. For example, 
IChabrier & Bar affe (1991) predict the boundary to occur at M « 0.35 M© (spectral type 



M2), whereas iDorman et al., (1989) place the fully convective boundary at M v 0.25 



Mq (spectral type M 4). For the fiducial propert ies of M2 stars, refer to the benchmark 
systems Gliese 22A (iHenrv & McCarth\l ll993') and Gliese 41 1 (iLane et al.ll200ll') . both 



with empirical masses of wO.36 M© and the latter with an interferometrically measured 
radius. For the fiducial properties of M4 stars, refer to the benchma rk echpsing binar y 
CM Dra, comprising two M4.5 stars with masses of 0.22 Mo (e.g. Moral es et all2009h . 
Table [Ugives a summary of the physical properties of spectral type M2 and M4.5 stars. 



Table 1 . Basic physical properties of exemplar M2 and M4.5 dwarfs. 



Name Mass (Mq) Radius (Rq) Teff My V-K 



GJ22A 0.352 + 0.036 ... 3380+ 150 10.99 + 0.06 4.55 

GJ411 0.403 + 0.020 0.393 + 0.008 3828 + 100 8.08 + 0.09 4.25 

CM Dra A 0.2310 + 0.0009 0.2534 + 0.0019 3130 + 70 12.78 + 0.05 5.10 

CM Dra B 0.2141 +0.0010 0.2396 + 0.0015 3120 + 70 12.91 +0.05 5.10 
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2.1. Basic stellar properties 

Neill Reid discussed the observed properties of M dwarfs across the fully convective 
boundary. It is useful to recall up front that certain fundamental stellar properties do not 
evince any obvious manifestation of the changes in stellar structure that accompany the 
transition into the fully convective regime. For example, the mass-luminosity relation is 
smooth through the M spectral types. Similarly, Figure [T^ shows the radius- luminosity 



relation for a large sample of M dwarfs (.Ribas .20061 : Demory et al.ll2009h : The ob- 
served radius-luminosity relation is smooth through the M spectral types, as predicted 
by stellar evolution models. Evidently, energy generation and output in these stars does 
not much know or care about changes in structure or in energy transport that may be 
occurring across the fully convective boundary. 



o.s . I . — ^ — ^ — ^ — ^ — t — ^ — ^ — t — t — t — t — ^ — t — t — t — ^ — ^ — ^ — ^ ^ 1 — r-i O.B 




Figure 1. Left: Radius-luminosity relation for M dwarfs. Crosses are nearby field 
M dwarfs with reliable parallaxes, and temperatures derived from V - I colors. The 
relationship is smooth across the M spectral types, as predicted by s tellar evolution 
models (dashed lines). Right: Radius-Teff relation. As highlighted ICle mens et al.l 
0998), there is a marked decrease in radius near spectral type M4. Note, however, 
that this change in the radius-Teff relationship is not infinitely sharp, but rather is a 
finite transition from spectral types M3 to M6. 

However, other basic stellar properties do manifest clear changes across a mass of 
~0.3 Mq. Figure [T}) shows the radius-temperature relation for the same stars in Fig- 
ure [T^. Unlike the featureless radius-luminosity relation, here we see a clear change in 
slope — with the stellar radii decreasing rapidly with decreasing effective temperature — 
centered around Teff ~ 3150 K (spectral type M4). The change in stellar radius here 
is a factor of ~2. In other words, there is a relatively large range of stellar mass and 
radius over a relatively small range of spectral types. Indeed, stars with a wide range 
of masses around the fully convective trans ition are predicted to have similar spectral 
types (cf. Fig. 4 in lChabrier & Baraffell2000l') . a consequence of the flatness of the mass- 
temperature relation caused by H2 formation. There is a strong morphological similar- 
ity between the observed and predicted relations, but the two are offset by ~250K, with 
the observations cooler than the models (compare pink and black curves in Fig.[T]). 

Notably, the change in the radius-temperature relation near M4 is not an abrupt 
step function; rather, it is a. finite transition between Tf-n * 3300 K and Teff « 2850 K 
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(i.e. spectral types M3 to M6). At later spectral types, the stellar radii asymptotically 
approach the value of wO. 1 Rq dictated by electron degneracy pressure. 

As Leslie Hebb discussed, recent direct radius measurements of M dwarfs in 
eclipsing binaries show the stellar radii t o be inflated relative to the radii predicted by 
theoretical stellar evolution models (e.g. lRibail2006). Thi s has been attributed to the 
eff'ects of stellar magnetic fields and activity (e.g. [Morales et al. 20 08). The magnitude 
of this eff'ect does change across the fully convective boundary (Figure |2]), but as above 
it does not appear to change in an abrupt fashion. 




Figure 2. Left: Mass-radius relation for eclipsing binary systems with masses and 
radii measured to better tha n 3%. Fract ional increase in the radii over the 

theoretically predicted radii (Baraf fe et al.lll998h . The enlarged radii are thought to 
be caused by the affects of magnetic fields. Early type M dwarfs show a larger 
increase in radii compared to late-type, fully convective M dwarfs, consistent wit h 
predictions of theoretical models (.Mullan & MacDonaldl2001l;IChabrier et alj2007h . 



2.2. Stellar activity and stellar rotation: Effects of stellar age 

As Jonathan Irwin and Richard Jackson discussed, recent observational results have 
confirmed that the M4 transition plays an important part in shaping the rotation and 
magnetic activity evolution of M dwarfs. Large samples of M dwarfs show an increase 
in activity fractions around spectral types of ~M4 dWest et al...2008.) . 




MO Ml M2 M3 M4 M5 UB MT 
Spectra! Type 



Figure 3. Activity lifetime as a function of spectral type for M dwarfs (West et 
al. 2008). The lifetime over which an M dwarf stays active is a strong function of 
spectral type, from 1-2 Gyr for M3 dwarfs to 7-8 Gyr for M6 dwarfs. Evidently, 
activity lifetime does respond to the change in internal structure across the fully 
convective boundary, but the transition is a finite one between M3 and M6. 
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This result has been shown to be dependent on the ages of the stars. When looking 
at the activity fractions as a function of the M dwarfs' location above the Galactic 
plane, and using a simple model of dynamica l heatin g which scatters M dwarfs to larger 
Galactic scale heights over time, IWest et al.l (l2008l) determine the amount of time that 
M dwarfs of different spectral types remain chromospherically active. They find that 
the "activity lifetime" changes dramatically between spectral types M3 and M6 (from 
1-2 Gyr for M3 dwarfs to 7-8 Gyr for M6 dwarfs; FigEll. 

This important transition was bolstered by the rotation analysis of lReiners & Basri 

(l2008h . who found that the average rotation velocities of a nearby sample of M dwarfs 
increases dramatically and monotonically as one looks at later spectral types. The av- 
erage vsin/ is <3 km/s at M2, increases to ~10 km/s at M6, and continues increasing 
to ~30 km/s at L4 (Figure |4l). This almost certainly points to a lengthening of the 
timescale for stellar spindown as the stellar structure changes from solar-type (i.e. ra- 
diative core -I- convective envelope) to fully convective. However, the transition is not a 
step function, thus again pointing to the importance of evolution over time. 




MO M2 M4 M6 M8 LO L2 L4 L6 L8 TO T2 T4 T6 T8 
Spectral Type 



Figure 4. Rotation (v sin /) of late-type stars (iReiners & Basrill2008h . Starting at 
spectral type ~M3, mean v sin / increases monotonically, implying that the spindown 
timescale increases monotonically across the fully convective boundary. 



The importance of time evolution is seen most clearly in rotation-period distri- 
butions in populations of different ages. Fig. [5] shows the evolution of stellar rotation 
periods from ~1 Myr to >10 Gyr. The situation among the pre-main-sequen ce popu- 
lations is complex, likely due to the effects of multiple competing phenomena (IBouvier 
120091 However, at ages of ~100 Myr and older, the stars show a clear pattern with 
two dominant sequences of rapid and slow rotators; the slow-rotator sequence becomes 
increasingly dominant with increasing stellar age. T hese are the so-calle d 'C and T 
sequences of the Barnes "gyrochronology" paradigm (lBamesll2003l . 120071) . 

Importantly, it is also clear from this figure that at all ages there is a spread of 
stellar rotation periods near the fully convective boundary; age alone does not determine 
whether a given M4 star will be on the 'C' sequence or on the T sequence. Thus, while 
all stars appear to follow a similar pattern of evolution in stellar rotation, presumably 
connected to evolution of their interior structure and magnetic field properties, there is 
evidently star-to-star scatter in these properties at any given age. 



2.3. Magnetic fields: Topology and effects on stellar properties 

As discussed by Moira Jardine and Julien Morin, M dwarfs also show clear topological 
changes in their surface magnetic fields as a function of spectral type. Figure [6^ shows 
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Figure 5. Observed rotation period distributions for low-mass stars in clusters of 
varying ages. Among clusters with ages > 100 Myr, two sequences of rapid and slow 
rotators are observed, with the slow -rotator branch becoming increasingly dominant 
with time. At the fully convective boundary (vertical line), the stars do not transition 
to the slow-rotator sequence instantaneously, but rather on a finite timescale. 



a simple cartoon illustration summarizing the basic change that is observed: from glob- 
ally weak but highly structured (i.e. multipolar and non-axisymmetric) fields at early 
spectral types, to globally strong and well ordered (i.e. dipolar) fields at later spectral 
types. In other words, the energy in the dipole component, relative to higher-order 
components, is higher in the lower-mass stars. Figure[6j5 depicts these changes in refer- 
ence to three specific low-mass stars (Gl 4 94, EV Lac, EQ Peg b) whose magnetic field 
topologies have been modeled in detail dJardine et alJ 12002) iising Zeeman Dop pler 
Imaging (ZDI) maps as boundary conditions ("Ponati et al. 20081 : iMorin et alj|2008l) . 

This is shown more quantitatively in Figure|7J which depicts the field topology and 
strength as a function of stellar mass, derived for a large number of low mass stars using 
the ZDI technique. Interestingly, we see a hint of two transitions here. The first occurs 
near the fully convective boundary, where as already mentioned the fields transition 
from a globally weak, highly structured topology, to a globally strong, dipolar topol- 
ogy. A second transition appears at spectral type M6-M8, such that below M~ 0. 15 Mq, 
stars with similar stellar masses and rotational properties can have very different mag- 
netic topologies. This apparent M6-M8 transition could b e reflecting the point a t 
which the stellar atmospheres become completely neutral (e.g. Moh anty & Basrill2003 ). 
though it could again reflect an age effect as witnessed in the bifurcated rotational prop- 
erties of M dwarfs (see Fig.[5l). 

Finally, it is important to bear in mind that the ZDI measurements, because they 
rely upon observations of polarized light, are probing field components that represent 
only a few percent of the total magnetic energy. There is likely a large amount of mag- 
netic flux in small-scale field components (which do not contribute much net polarized 
light) in addition to the larger scale components that the ZDI technique is mapping. 
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Figure 6. Left: Cartoon representing observed differences in surface magnetic 
field topology across the fully convective regime. Right: Reconstructed surface field 
topologies for three low-mass stars are shown (top to bottom: Gl 494, EV Lac, EQ 
Peg b). The field changes from a strong, highly structured, multi-polar configuration 
at early M types, to a weaker, more well organized dipole at the late M types. 

Dermott MuUan discussed recent modeling efforts to account for the observed in- 
flated radii of M dwarfs compared to the predictions of most (non-magnetic) theoretical 
stellar evolution models. Chabrier and collaborators have suggested that strong surface 
fields can inhibit surface convection and also create cool starspots, both of which act to 
decrease the surface temperature and then — because the stellar luminosity generated in 
the core is unaffected — requires an increase in the stellar radius. 

Mullan also discussed alternate models in which a strong field threads the entire 
star. These models can successfully reproduce for example the observed temperature 
reversal in the brown-dwarf eclipsing binary 2M0535-05 (FigureO, which at an age of 
only ~ 1 Myr has a spectral type of M6. These models do predict a change in the stellar 
luminosity as a result of the very strong (~MG) fields in the core (Figure |9ll. 



3. Putting it all together: A conceptual framework 

Figure \T0\ is an initial attempt to conceptually summarize the inter-related phenomena 
discussed during our Session. This is not intended to be inclusive of all mechanisms or 
observable phenomena possibly relevant to M dwarfs, and it is by necessity a simplified 
representation of what is likely to be a more complex set of inter-relationships. 

The figure attempts to represent the complex set of causal relationships — stellar 
structure, magnetic field topology, rotation, magnetic field strength — that mediate the 
fundamental relationships between basic stellar properties (e.g., how mass and com- 
position ultimately translate into radius and effective temperature). As depicted in the 
figure, we imagine that, absent any magnetic field effects, there would be a single, fun- 
damental mass-radius-temperature relationship (modulo metallicity). This relationship 
is, however, modified by the presence of a magnetic field that manifests itself on the sur- 
face through various activity tracers (e.g. chromospheric Ha emission). This magnetic 
field affects the observed surface prop erties of the star (effe ctive temperature and ra- 
dius) via surface effects (i.e. star spots: IChabrier et al] (l2007t)) and/or via a strong field 
that inhibits convection throughout the stellar interior (Mullan & MacDonald 2001). 

The mechanism that drives the magnetic field's action (and perhaps its generation) 
is likely complex, but almost certainly involves at its heart a rotational dynamo. We 
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Figure 7. Properties of the large-scale magnetic topologies of a sample of M 
dwarfs inf erred from spectropolarimetric observations and Zeeman D oppler Imag- 
ing (ZDI) CiDonati et al.ll2008HMorin etani2008t|Phan-Bao et alj|2009 ) as a function 
of rotation period and mass. Larger symbols indicate stronger fields, symbol shapes 
depict the degree of axisymmetry of the reconstructed magnetic field (from decagons 
for purely axisymmetric to sharp stars for purely non axisymmetric), and colours the 
field configuration (from blue for purely toroidal to red for purely poloidal). Left- 
ward white arrows indicate stars for which only an upper limit on Pjot is derived and 
the field properties are extrapolated. For VB8 the same symbol as VBIO is arbi- 
trarily used, no magnetic field was detected from spectropolarimetric observations. 
For GJ 1245 B symbols corresponding to 2007 and 2008 data are superimposed in 
order to emphasize the variability of this object. Solid lines represent contours of 
constant Rossby number Rq = 0.1 (saturation threshold) and 0.01. The theoretical 
full-convection hmit (M~ 0.35 Mq, Chabrier & Baraffe ( 1997)) is plotted as a hor- 
izontal dashed line, and the approximate limits of the three stellar groups identified 
are represented as horizontal dotted lines. 



envision that the star's rotation directly affects the strength of the magnetic field (more 
rapid rotation — > stronger field) for stars below the saturation threshold. In general, a 
stronger field will drive a more rapid evolution of the star's angular momentum, leading 
to a decline of the star's rotation on a certain spindown timescale, and this spindown in 
turn feeds back into the strength of the field, dialing it down over time. 
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Figure 8. The reversal of effective temperature with mass for the brown-dwarf 
ecHpsing binary 2M0535-05 ( iStassun et alj|2006[ 2007) can be explained by recent 
models that include the effec ts of strong magnetic fields th reading through the entire 
interior of a low-mass star (iMacDonald & Mullanll2009h . Solid curves show the 
non-magnetic models for the primary (more massive) brown dwarf (red) and the 
secondary brown dwarf (blue), dashed lines show the magnetic models. 



Importantly, the magnetic field acts according to an intrictate feedback process 
such that the magnetic field impacts the stellar structure (i.e. convection zone depth) and 
rotation rate of the star which in turn affects the (re)generation of the field. The nature 
of the dynamo itself depends on the interior structure (i.e., convection zone depth). As 
the star approaches a more fully convective state (that is, as the size of the radiative 
zone recedes), there must be a transition from a solar- type (a - oj) dominated dynamo 
to a purely turbulent (a^) dynamo. The spindown timescale is probably intimately 
coupled to the field topology, but exactly how the two are connected is not obvious. 
For example, later type M dwarfs have longer spindown timescales, yet these same 
stars have magnetic field topologies which posses more open field lines (i.e. dipolar 
fields), which would seem to enable them to more effectively drive mass and angular 
momentum loss. Clearly, we do not yet understand entirely how magnetic topology 
affects the observed rotational evolution of these stars. 

It is moreover not clear whether the dynamo-driven strength of the field can feed 
back into the depth of the convection zone and/or the field topology. For example, it is 
thought that a very strong field can cause an otherwise fully convective star to generate 
a radiativ e core, either directly throug h the inhibition of convection at the center of the 



star (e.g. MuUan & MacDo naldl 12001) and/or indirectly th rough altering the boundary 



conditions at the stellar surface (e.g. IChabrier et aLll2007h . Similarly, in some models 
very r apid stellar rotation can cause the opening of surface field lines (e.g. Ijardine et al.l 



2010f) . Thus, the inter-relationships between rotation, field strength, interior structure. 



and surface field topology are likely complex, non-linear, and in any event dependent 
on the age of the star through its rotational history. 

The rotation-period distributions of low-mass stars as a function of stellar age show 
that rotation is not simply a function of time alone. Indeed, current thinking identifies 
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Figure 9. Models of iMullan & MacDonaldl ( 1200 Ih showing the change in stellar 
radius, temperature, and luminosity, resulting from strong internal magnetic fields. 



at least two dominant sequences of rotational properties in clusters (e.g., the 'C and T 
sequences of Barnes), such that a star of a given mass and age is likely to have one of 
two possible rotation periods. Perhaps these distinct rotational sequences reflect stars in 
different states of the complex feedback mechanisms between field generation, dynamo 
type, and field topology explored above. Indeed, the ZDI maps of field topologies ap- 
pear to bear out this picture of multiple types of field configurations intermixed among 
the M3-M6 spectral types (Figure |7]). Given this, and in light of the time-dependence 
of these effects, it is perhaps not surprising that the "M4 transition" across the fully 
convective boundary is not a sharp transition, but rather "smeared" from approximately 
M3 to M6. For example, the two eclipsing binaries CM Dra and CU Cnc (see Fig. [T|l, 
with masses and radii that differ by a factor of ~2 but with nearly identical effective 
temperatures (spectral types wM4), may be exemplars of this smearing effect. 

4. Open questions and future work 

We concluded with an open discussion leading to a set of questions to help guide ongo- 
ing work by the community interested in the "M4 transition" and related phenomena. 

1. What field strengths are necessary and reasonable in the interiors of stars to po- 
tentially explain the observed mass-radius relationship? 

2. Abundance measurements needed! (e.g. [0/Fe]) The mass-radius relationship is 
actually a mass-radius-abundance relationship. FGK-i-M binaries are a promising 
approach, tying M dwarf abundances to FGK primary stars. 
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Figure 10. Representation of the process by which magnetic fields influence the 
interior structure and observed fundamental properties of M dwarf stars. M dwarfs 
are defined intrinsically by mass and chemical composition which determine the ini- 
tial depth of the convection zone. This in turn affects the magnetic field topology and 
strength through the dynamo type (a-Q. or a^). The rotation rate, which is moder- 
ated by angular momentum loss through stellar winds determined largely by the field 
topology, in turn affects the field strength for stars in the non-saturated regime via the 
dynamo process. X-rays, chromospheric emission, starspots, and flaring all trace the 
magnetic activity. Over time, stars spin down and the magnetic field weakens, caus- 
ing the entire system to evolve. It remains unclear how much feedback there is in the 
system. Does the strength of the magnetic field affect the topology or the convection 
zone depth, and how much does the rotation rate influence the field topology? Taken 
together, the interlinked magnetic fields, internal structure, and rotation properties 
affect the resulting fundamental parameters: radius and temperature. However, the 
detailed mechanism by which the latter occurs is not yet fully understood. 



3. Continue to identify low-mass eclipsing binaries. Long-period systems are needed 
to disentangle effects of magnetic activity on stellar radii and temperatures. 

4. Examine selection effects in samples, especially with regard to activity. Ages of 
stellar samples matters. 

5. What is the mass of the fully convective boundary? Depends on rotation/activity/age? 

6. Other observational clues: Period gap in CVs (2-3 hr periods) — > kink in mass- 
radius relationship? 

7. Independent radius measurements: interferometry, P^t + v sin / (gives R sin /) 
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8. Bulk of field is in very small-scale structures — > How to fully characterize the 
fields? Constraints from many techniques/tracers are needed to fully characterize 
magnetic fields: Zeeman broadening, spectropolarimetry, radio, X-rays, etc. 

9. "M6-M8 transition" in field topologies, caused by neutral atmospheres? Age? 

10. How does the fully- to partly-convective evolution on the pre-main-sequence 
affect the early magnetic and rotational history of low-mass stars? 
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